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ABSTRACT 

The mechanism by which outflows and plausible jets are driven from black hole systems, still re- 
mains observationally elusive. Notwithstanding, several observational evidences and deeper theoretical 
insights reveal that accretion and outflow/jet are strongly correlated. We model an advective disk- 
outflow coupled dynamics, incorporating explicitly the vertical flux. Inter-connecting dynamics of 
outflow and accretion essentially upholds the conservation laws. We investigate the properties of the 
disk-outflow surface and its strong dependence on the rotation parameter of the black hole. The 
energetics of the disk-outflow strongly depend on the mass, accretion rate and spin of the black holes. 
The model clearly shows that the outflow power extracted from the disk increases strongly with the 
spin of the black hole, inferring that the power of the observed astrophysical jets has a proportional 
correspondence with the spin of the central object. In case of blazars (BL Lacs and Flat Spectrum 
Radio Quasars), most of their emission are believed to be originated from their jets. It is observed 
that BL Lacs are relatively low luminous than Flat Spectrum Radio Quasars (FSRQs). The lumi- 
nosity might be linked to the power of the jet, which in turn reflects that the nuclear regions of the 
BL Lac objects have a relatively low spinning black hole compared to that in the case of FSRQ. If 
the extreme gravity is the source to power strong outflows and jets, then the spin of the black hole, 
perhaps, might be the fundamental parameter to account for the observed astrophysical processes in 
an accretion powered system. 

Subject headings: accretion, accretion disks — black hole physics — hydrodynamics — galaxies: active 
— galaxies: jets — X-rays: binaries 



1. INTRODUCTION 

High resolution observations show strong outflows and 
jets in black hole accreting systems, both in active galac- 
tic nuclei (AGNs) or quasars (Begelman et al. 1984; 
Mirabel 2003) and microquasars (SS433, GRS 1915+105) 
[Margon 1984; Mirabel & Rodriguez 1994, 1998]. Ex- 
tragalactic radio sources show evidence of strong jets in 
the vicinity of spinning black holes (Meier et al. 2001; 
Meier 2002). Outflows are also observed in neutron star 
low mass X-ray binaries (LMXBs) (Fender et al. 2004; 
Migliari & Fender, 2006) and also in young stellar objects 
(Mundt 1985). It has been argued (Ghosh & Mukhopad- 
hyay 2009; Ghosh et al. 2010; hereinafter GM09, G10 
respectively, and references therein) that outflows and 
jets are more prone to emanate from strong advective 
accretion flows; the said paradigm is more susceptible 
for super-Eddington and sub-Eddington accretion flows. 
However, the exact (global if any) mechanism of forma- 
tion of the jet, its collimation, acceleration, composition 
from the accretion powered systems still remain incon- 
clusive. 

Extensive works have been pursued on the origin of 
outflow/jet, since the pioneering work of Blandford & 
Payne (1982) (e.g. Pudritz & Norman 1986; Contopou- 
los 1995; Ostriker 1997), where the authors used a self- 
similar approach to demonstrate that the poloidal com- 
ponent of the magnetic field can be seemingly used to 
launch outflowing matter from the disk. The formation 
of the jet is directly related to the efficacy of extrac- 
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tion of angular momentum and energy from the accre- 
tion disk. Physical understanding of the hydromagnetic 
outflows from disks has been developed from magnetohy- 
drodynamic (MHD) simulations (Shibata & Uchida 1986; 
Ustyogova et al. 1999; De Villiers et al. 2005; Haw- 
ley & Krolik 2006) both in non-relativistic as well as in 
relativistic regimes, mostly in the Keplerian paradigm. 
However, the large Lorentz factors as well as Faraday 
rotation measures suggest that the observed VLBI jets 
in quasars and active galaxies are in the Poynting flux 
regime (Homan et al. 2001; Lyutikov 2006). On the 
other front, strong radiation pressure can serve as a dif- 
ferent mechanism to effuse outflows /jets. This is likely 
to occur when the accretion rate is super-Eddington 
or super-critical (Fabrika 2004; Begelman et al. 2006; 
GM09,G10 and references therein) and the accretion disk 
is precisely "radiation trapped" as in ultra-luminous X- 
ray (ULX) sources. It has been confirmed by radiation 
hydrodynamic simulation at super-critical accretion rate 
(Okuda et al. 2009) as well to explain luminosity and 
mass outflow rate of relativistic outflows from SS433. 
The under-luminous accreting sources, having high sub- 
critical accretion flow, were explained by an advection 
dominated accretion flow (ADAF) model (Narayan & Yi 
1994). The promising outcome of this model lies in the 
large positive value of the Bernoulli parameter because 
of the small radiative energy loss. It leads to conceive 
that the gas in the inflowing disk is susceptible to escape, 
leading to strong unbounded flows in the form of outflows 
and jets. This also signifies that even in the absence of 
magnetic field and radiation pressure, outflows are plau- 
sible from strongly advective accretion flow if the system 
is allowed to perturb. Nevertheless, the definitive under- 
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standing of the origin of outflows/jets is sill unknown. 
It remains one of the most compelling problems in high 
energy astrophysics. 

Whatever might be the reason for the origin of out- 
flow and then jet from the disk, one aspect is however 
definite that strong outflows producing relativistic jets 
are powered by extreme gravity. Although seems para- 
doxical, the strength and the length scale of observed 
astrophysical jets vary directly with the strength of the 
central gravitating potential. Observationally, it is evi- 
dent that strong outflows and relativistic jets are more 
powerful in observed AGNs and quasars, harboring su- 
pcrmassive black holes, compared to that in black hole 
X-ray binaries (XRBs). In addition, the jets observed in 
AGNs and quasars have greater length scale compared to 
that seen in stellar mass black hole systems. One of the 
most important signatures of relativistic gravitation is 
the spin of the central object. It is presumably believed 
that the spin (practically specific angular momentum) of 
the neutron star is less than that of a black hole, and 
thus the observed jet from black holes is much stronger 
and powerful than that of neutron star sources. In early, 
Blandford & Znajek (1977) demonstrated that if there 
is a magnetic field associated with the black hole due 
to threading of magnetic field lines from the disk and 
the angular momentum of the Kerr black hole is large 
enough, then the energy and the angular momentum can 
be extracted from the underlying black hole by a purely 
electromagnetic mechanism, which can thus be expected 
to power the jet in an AGN. These imply that the spin 
might play a significant role in powering jets, both in 
microquasars and in AGNs, rather, it can act as a fun- 
damental parameter in an accretion powered system. 

Most of the studies of accretion disk and studies of re- 
lated outflow /jet have been evolved separately, assuming 
these two apparently to be dissimilar objects. However, 
several observational inferences (for details see GM09, 
G10 and references therein) and improved understand- 
ing of accretion flow and outflow reveal that accretion 
and outflow/jet are strongly correlated. The unifying 
scheme of disk and outflow is essentially governed by con- 
servation laws; conservations of matter, energy and mo- 
mentum. Hence, in modeling the accretion and outflow 
simultaneously in any accretion powered system, follow- 
ing aspects should be taken into account: (1) the effect 
of relativistic central gravitational potential including its 
spin, (2) proper mechanism of origin of outflow/jet from 
the disk, (3) appropriate hydrodynamic equations (con- 
sidering that the accreting gas be treated as a continuum 
fluid), capturing the information about the intrinsic cou- 
pling between inflow and outflow which are governed by 
the conservation laws in a strongly advective paradigm. 

Recently, GM09, G10 made an endeavor to explore the 
dynamics of the accretion-induced outflow around black 
holes/compact objects in a 2.5-dimensional paradigm. 
The authors formulated the disk-outflow coupled model 
in a more self-consistent way by solving a complete set of 
coupled partial differential hydrodynamic equations in a 
general advective regime through a self-similar approach 
in an axisymmetric, cylindrical coordinate system. They 
explicitly incorporated the information of the vertical 
flux in their model. However, they restricted their study 
to Newtonian regime, thus negating the requisite effect of 
general-relativity, especially the effect of spin of the cen- 



tral object. Based on the model, the authors computed 
the mass outflow rate and the power extracted by the 
outflow from the disk self-consistently, without propos- 
ing any prior relation between the inflow and outflow. 

In the present paper, we propose a new model for 
the accretion-induced outflow by extending the work of 
GM09, G10, by incorporating the general relativistic ef- 
fect of the central potential without limiting ourselves 
to a self-similar regime. As the definitive mechanism of 
launching of outflows/jets is still evasive, we do not em- 
brace any specific mechanism of outflow like inclusion of 
the magnetic field into our model equations. Neverthe- 
less, the importance of the magnetic field can not be, 
in principle, discarded to explain the launching and col- 
limation of jets off the accretion disk. Here, owing to 
our inability to solve coupled partial differential MHD 
equations in a 2.5-dimensional advective regime, we ne- 
glect the influence of the magnetic field in our model 
equations. However, the implicit coupling between the 
inflow and outflow, dictated by the conservation equa- 
tions, have been taken into account appropriately. We 
arrange our paper as follows. In the next section, we 
present the formulation of our model. The section 3 de- 
scribes the computational procedure to solve the model 
equations of the accretion-induced outflow. In sections 4 
and 5, we study the dynamics and the energetics of the 
flow respectively. Finally, we end up in section 6 with a 
summary and discussion. 

2. MODELING THE CORRELATED DISK-OUTFLOW 
SYSTEM 

We formulate the disk-outflow coupled model by con- 
sidering a geometrically thick accretion disk, which is 
strongly advective as strong outflows /jets are more likely 
to eject from a thick/puffed up region of the accretion 
flow (GM09,G10). The vertical flow is explicitly included 
in the system. The basic features of the model are similar 
to that in GM09. We adopt the cylindrical coordinate 
system to describe a steady, axisymmetric accretion flow. 
The dynamical flow parameters, namely, radial velocity 
(v r ), specific angular momentum (A), vertical velocity or 
outflow velocity (v z ), adiabatic sound speed (c s ), mass 
density {p) and pressure (P) depend both on radial and 
vertical coordinates. We have already highlighted the 
importance of the spin of black hole to power the out- 
flow/jet and its presumed relative effect on the obser- 
vation of various AGN classes. Thus we have included 
the effect of the spin in our model. As the spin of the 
black hole is a signature of pure general relativity, i.e., 
Einstein's gravitation, its effect on the accretion flow, 
especially in the inner region of the disk, is mimicked ap- 
proximately with the use of pseudo-general-relativistic 
or pseudo- Newtonian potential (PNP). Because of the 
disk-outflow system to be geometrically thick and the 
flow to be 2.5-dimensional (not confined to the equato- 
rial plane), we use the PNP of Ghosh & Mukhopadhyay 
(2007), which is a pseudo-Newtonian vector potential, to 
capture the inner disk properties of the accretion flow 
around a Kerr black hole approximately. 

At the first instant we neglect the effect of viscosity in 
our system. One of the reasons behind it is the unavail- 
ability of the effective computational technique to solve 
coupled partial differential viscous hydrodynamic conser- 
vation equations for the compressible flow. To make the 
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inviscid assumption more arguable, it can also be noted 
that the angular momentum transport in the accretion 
flow, for which the necessity of turbulent viscosity is in- 
voked, can alternatively take place purely by outflow. At 
this extreme end, the outflow extracts angular momen- 
tum from the disk allowing the matter to get accreted 
towards the black hole and hence the inviscid assump- 
tion can be adopted. This is in essence similar to the 
Blandford & Payne (1982) mechanism to extract energy 
and angular momentum from the magnetized disk, where 
the extraction of the angular momentum and energy is 
essentially done by the outflow, and not due to the vis- 
cous dissipation. The outflow originates from just above 
the equatorial plane of the disk and this lower boundary 
is maintained at z = 0, when v z = 0, unlike other work 
where the outflow is hypothesized to effuse out from the 
disk surface (e.g. Xie & Yuan 2008). Our model is ef- 
fectively valid only in the region where disk and outflow 
are coupled, i.e. the region from where essentially the 
outflow is emanated from the accretion flow. Hence our 
study will remain confined within this predefined region. 
Further, we neglect the contribution of the magnetic field 
as argued in §1. 

We circumvent the idea of vertical integration of the 
flow equations. The validity and the reliability of the 
height integrated equations is normally gratifying in the 
geometrically thin limit. In that circumstances, the flow 
velocities are likely to be more or less independent of the 
disk scale-height, which is not the case of the present 
paradigm of interest. We further consider that the disk 
to be non-self-gravitating, assuming that the mass of the 
disk to be much less than that of the black hole. The 
radial and vertical coordinates are expressed in units of 
GM/c 2 , flow velocities in c, time in GM/c 3 and the spe- 
cific angular momentum in the unit of GM/c. Here G, 
M and c are gravitational constant, mass of the black 
hole and speed of light respectively. The steady state, 
axisymmetric disk-outflow coupled equations in cylindri- 
cal geometry in the inviscid limit are then as follows. 

(a) Mass transfer: 

Id, . d , . , . 

-^-{rpvr) + -^-{pvz) = 0, (1) 
r or oz 

where the first term is the signature of accretion and the 

second term attributes to outflow. As the outflow starts 

from just above z — surface, within the inflow region 

itself, we make a reasonable hypothesis that within the 

prescribed disk region the variation of the dynamical flow 

parameters with z is much less than that with r, allowing 

us to choose dA/dz w 0(A/z); for any parameter A. 

Strictly speaking, the weak variation with z ensures that 

the outflow originating from the mid-plane of the disk 

does not disrupt the disk structure, and thus allowing a 

smooth accretion flow towards the black hole. Thus Eqn. 

|T]) then reduces to 

' % Ur ) + ^ = 0. (2) 



r dr 

(b) Radial momentum balance: 

dv r dv r A 2 

v r ~z 1- v z — + t Gr 

Or oz r^ 



if = 0, (3) 
p or 

where Fq t is radial component of the gravitational force. 
As discussed earlier, here Fc r is the radial force corre- 
sponding to the PNP in cylindrical coordinate system 



given by Ghosh & Mukhopadhyay (2007) containing the 
information of spin of the black hole. With the similar 
argument as above the equation then reduces to 



8v r V r -V r o A 2 1 BP 

v r — h v z T + F Gr H — = 0, (4) 

z r° p Or 



dr 



where v r g is the radial velocity at the mid-plane of the 
disk. 

(c) Azimuthal momentum balance: 



1 d i 2 s 



d 

-q- z (pv<pv z ) = 0, 



(5) 



where v$ is the azimuthal velocity of the flow. The first 
term of this equation signifies the radial transport of the 
angular momentum, while the second term describes the 
extraction of angular momentum due to mass loss in the 
outflow. If we consider that the net angular momentum 
extracted by the outflow be Xj , and the remaining angu- 
lar momentum retained by the disk A<2, then total angu- 
lar momentum A = \ 3 ■ + A^ can be assumed to remain 
constant throughout the flow within our predefined disk- 
outflow region by the virtue of an inviscid flow. There- 
fore, 



A = constant, 
(d) Vertical momentum balance: 



dv z 
dr 



dv z 
dz 



F G z 



1 dP 

p dz 



0, 



(6) 



(7) 



where Fq z is the vertical component of the gravitational 
force, described by Ghosh & Mukhopadhyay (2007). Fol- 
lowing previous arguments this reduces to 



dv z 
dr 



z 



F 



Gz 



l P-Po 

p z 



= 0, 



(8) 



where Pq is the pressure of the flow at the mid-plane of 
the disk. If there is no outflow, then v z = 0, and Eqn. © 
reduces to the well known hydrostatic equilibrium condi- 
tion in the disk, and the hydrostatic disk scale-height can 
be obtained. Similarly, one can customarily extend the 
vertical momentum equation to compute the disk scale- 
height when there is an outflow coupled with the disk. 
Thus we will use Eqn. (JSJ) to obtain the scale-height 
of the disk-outflow coupled system. We can reasonably 
assume that at height h (i.e. at z—h), the pressure of 
the disk is much less compared to that at the equatorial 
plane to prevent any disruption of the disk, permitting 
a steady structure of accretion flow. The variation of 
density along z direction is approximately kept constant 
(po ~ p) which in turn means that the disk-outflow cou- 
pled region is weakly stratified. Considering the above 
facts, Eqn. © is simplified to obtain the scale-height as 



i dv z 



Vl\h 



F, 



Gz\h 



hp 



(9) 



where h is thus the the root of the above transcendental 
equation. Thus Eqns. will simultane- 

ously have to be solved to obtain the dynamics and the 
energetics of the accretion-induced outflow. 
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3. SOLUTION PROCEDURE AND THE DISK-OUTFLOW 
SURFACE 

We assume that the accretion flow follows an adiabatic 
equation of state P = Kp 1+1 / n where n — 1/(7 — 1), n 
and 7 are the polytropic and the adiabatic indices re- 
spectively. Note that constant K carries the informa- 
tion of entropy (e.g. Mukhopadhyay 2003) of the flow. 
Thus for an adiabatic flow, c s — (7^) 1 ^ 2 - In an ear- 
lier work (GM09), a reasonable assumption was made 
that h ~ r/2, which can be approximately accepted for a 
geometrically thick, 2.5-dimensional disk structure. Pre- 
suming that the outflow velocity is not likely to exceed 
the sound speed at the disk-outflow surface (the outer 
boundary of the accretion flow in the z direction), we 
propose a simplified relation between c s and v z as 



v z <2(-)c s 
r 



(10) 



Hence, at z = h ~ r/2, we obtain 



v z <c s . (11) 

This implies that the outflow can ideally come out off the 
disk surface which can appropriately be termed as the 
sonic surface in the vertical direction. With this notion 
in mind let us generalize this particular scaling between 
v z and c s , pertaining to our formalism as 



v z = i(~Tc s 
r 



(12) 



where 1 and /z are the constants which will be deter- 
mined by substituting Eqn. f| 12[) in our model conser- 
vation equations described in §2. The index \x measures 
the degree of subsonic nature of the vertical flow within 
the disk-outflow coupled region (i.e. in between the mid- 
plane and the surface of the accretion flow) . 

Generalizing the procedure adopted in earlier works 
for 1.5-dimensional disks (Chakrabarti 1996; Mukhopad- 
hyay 2003; Mukhopadhyay & Ghosh 2003), we solve 
Eqns. ©, (gj, © & dSJ. Using Eqn. jUJ along with 
Eqn. ((6J and combining Eqns. ^ & (j3J), we obtain 

dv r _ £-F Gr + & + I (f YCsjVrQ - V T + j|) _ N 

(13) 

Equation (fT3| shows that to guarantee a smooth solution 
at the "critical point", N = D = 0. From the critical 
point condition we obtain v rc — c sc , at r — r c . Here 
subscript c is referred to critical point. The radius r c is 
also called the "sonic radius" since no disturbance cre- 
ated within this radius can cross the radius (also known 
as sound horizon) and escape to infinity. Conditions at 
critical/sonic radius give 



Iz r r 



— — c sQc r c 
2z \ r r 



-r c Fcr 



A 2 . 



21 1/2 



(14) 



where c s0c = v r0c = (r c F Gr0c - A 2 / V 2 ) 1 ' 2 , is the sound 
speed or the radial velocity of the flow at the sonic ra- 
dius in the disk mid-plane (z = 0). F Grc is the radial 



gravitational force at the sonic radius and FcrOc is the 
corresponding value in the disk mid-plane. 

Now at sonic location = 0/0. Hence by applying 
['Hospital's rule, Eqn. (fT3|) reduces to 



dv r 
dr 



2n 



-B + y/B 2 -AAC \ vrc 
2A ) + r c 



+ - 



(15) 



where A = J"i(r, z)\ c , B = F%{r, z)\ c and C = J" 3 (r, z)\ c , 
are complicated functions of r and z at sonic location; 
J- 1, J- 2 & ^3 have been explicitly given in the appendix. 
Equations ([r3]) and ([15]) are then solved with an appro- 
priate boundary condition to obtain v r and c s as func- 
tions of r. The value of specific angular momentum in 
our flow always remains constant and is same as A c , the 
value at sonic radius, by virtue of Eqn. ©. 

Until now, we have emphasized on velocity profiles of 
the accretion-induced outflow at any arbitrary z. We 
have, however, mentioned before that the intrinsic cou- 
pling of inflow and outflow is confined within a speci- 
fied region, from where the outflow emanates. The disk- 
outflow inter-correlated region is bounded vertically from 
the mid-plane (z = 0) to an upper surface above which 
any inflow of matter ceased to exist. We aim at pre- 
cisely investigating the nature and dynamics of the flow 
at different layers within this disk-outflow coupled region. 
Therefore, we first need to calculate the scale-height of 
the accretion flow, based on our proposed model, and 
obtain the disk-outflow surface (upper boundary). It is 
to be noted that the transcendental Eqn. ([5]) cannot be 
independently used to calculate the scale-height h or its 
general variance with r, for unknown variables v r and c s . 
Nevertheless, we can easily obtain the scale-height h at 
r c , say h c , from Eqn. ([9]) as radial velocity and sound 
speed at the sonic location are known. Thus at r c , Eqn. 
((9|) simplifies to 



r c 



Csc\h £sc\h 

A* 



h c 

+Fc Z c\h 



dc s 
dr 



SOC Q 

"fh c 



(16) 



where dc s /dr\ c h is dc s /dr at the sonic location in the 
plane with scale-height h. Then h c obtained from Eqn. 
(fTB|) can be inserted into Eqn. (fTSj) to obtain dv r /dr\ c 
at scale- height h. 

The disk scale-height h without any outflow is seen to 
be linearly increasing with r (approximately by dimen- 
sional analysis). For a 2.5-dimensional disk-outflow sys- 
tem, it seems that h and r can also be linearly connected, 
based on the order of magnitude analysis (GM09). How- 
ever, unlike the thin disk, strong gas pressure gradient in 
the disk-outflow coupled system leads to a geometrically 
thick 2.5-dimensional disk. Therefore, for the present 
purpose we make a generalized scaling of h with r as 

h ~ 5r, (17) 

where S be a dimensionless arbitrary constant or any 
numerical variable (discussed in detail later). The im- 
pressionable choice about the factor 5 would be that it 
should contain precisely the information of the nature of 
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the flow. The only physical information we can extract 
related to the scale-height from our model equations is 
h c . Thus we rationally demand an approximate expres- 
sion for dimensionless parameter 5 as 

6 ~ ^. (18) 

r c 

It is found that the value of S circumvents around 1/2, 
for the entire range of spin parameter of the black hole, 
which here acts as a normalization constant. Thus even- 
tually, Eqns. (H3J), ([T5]). (IT71) & (0 are combined to- 
gether to obtain v r and c s along the scale-height h for an 
accretion-induced outflow. It also appears that under all 
these circumstances, and for a physically acceptable flow, 
the constant i and the index /i of Eqn. (|12p connecting 
v z and c s yield to be ~ 1 and ~ 3/2 respectively. 

3.1. Construction of disk- outflow surface 

After establishing the scale-height h, we in princi- 
ple can obtain the profiles of the dynamical parameters 
at different layers of the disk, i.e. at z = £h, where 
< £ < 1, within our prescribed disk-outflow region. We 
find that the radial velocity profile v r along all layers of 
disk for any arbitrary spin parameter exhibits some un- 
usual, yet very interesting behavior. For any specific z 
and spin parameter a, v r attains a negative value at a 
radius greater than a certain distance r — r\>. The mag- 
nitude of rf, decreases with the increase of £, which indi- 
cates that the positive trait of v r is greater for lower z. 
The negative value of v r does not necessarily reflect any 
unphysical behavior. If the outflow originates from any 
particular layer ih at a radius r > rt,, it then appears that 
the coupling between the disk and the outflow ceased to 
exist. It infers that along the layer £h, both accretion 
and outflow simultaneously occur up to the radial dis- 
tance ffe, and beyond which the solution (with negative 
v r ) reflects the truncation of the disk (along the specified 
layer). Let us consider different layers of z. As we as- 
cend to layers from a lower to a higher z, the truncation 
of the disk-outflow region along these layers occurs at a 
radius smaller than that of the lower z. As accretion is 
the source of the outflow, we can ostensibly conclude that 
the region from where the outflowing matter originates in 
the disk shrinks as we go to higher latitudes. Note that r& 
corresponds to zero v r . Therefore, we restrict our study 
up to the boundary r&. With this insight, we simply com- 
pute rf, along various layers of the disk and construct a 
surface in the r — z plane connecting all for different 
layers. The enclosed region bounded by this surface is 
defined by the positivity of v r , where both the accretion 
and outflow simultaneously persist and are intrinsically 
coupled to each other. We attribute the outer surface of 
the region as disk-outflow surface (h sur f) and above this 
layer no inflow takes place. The arrows in the diagram 
shown in Fig. 1 reveal the direction of flow. The arrow 
exactly at ri, points vertically upwards, which indicates 
that just at the flow pattern exclusively corresponds 
to vertical motion and any accretion flow ceased to ex- 
ist. However, more realistic disk-outflow surface could 
be visualized with a thick-solid line shown in Fig. 1. 

As in our system we have neglected the viscosity and 
any radiative loss, and thus the flow is considered to be 
strongly advective, presumably it is gas pressure dom- 
inated. Notwithstanding, various scattering processes 
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Fig. 1. — Nature and geometry of the disk-outflow coupled 
region and the outflow surface. Discussed in detail in §3.1. 
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Fig. 2. — Variation of disk-outflow surface with radial coordinate 
for different spin of the black hole. Solid, dashed, dot-dashed and 
dotted curves are for a = 0, 0.5, 0.9, 0.998 respectively. 7 = 1.5. 

will indeed produce radiation in the system, whose con- 
tribution could be insignificant. Hence, 7 ~ 3/2 is rea- 
sonably an appropriate choice in our model . Note that 
this is in essence similar to the choice of the earlier au- 
thors who modeled gas dominated low mass advection 
dominated accretion flows (e.g. Narayan & Yi 1994). 
We would also like to clarify that 7 = 5/3 corresponds 
only to the case of zero angular momentum. Figure 2 
shows the profiles of disk-outflow surface for various spin 
parameters a of the black hole. All the input parameters 
corresponding to each value of a are listed in Table 1 . It 
is seen that the disk-outflow region and the peak of the 
surface namely, Rj S , shifts to the vicinity of the black 
hole for a higher spin. 

The marginally stable orbit in a Keplerian accretion 
flow defines the inner edge of the disk (with zero torque 
condition). However the orbit is more or less an arti- 
fact of the flow. For a transonic advective disk (as is 

4 The relationship between 7 and /3, the ratio of gas pressure to 
the total pressure, is given by /3 = ^5§y (GM09). 
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Table 1 

Input parameters for various a 



a 


r c 


A 


0.0 


6.15 


3.3 


0.3 


5.7 


2.9 


0.5 


5.5 


2.6 


0.7 


5.3 


2.3 


0.9 


5.0 


1.9 


0.95 


4.4 


1.6 


0.97 


4.3 


1.5 


0.998 


4.2 


1.3 



the present case), an apriori definition of inner edge is 
somewhat fuzzy, due to which we can pretend to choose 
that the disk extends to the black hole horizon. In pres- 
ence of an outflow intrinsically coupled to the disk, there 
is supposed to be an inner boundary beyond which any 
outflow and then jet will ceased to exist. We define this 
inner boundary of the disk-outflow coupled region in ex- 
plaining the energetics of the flow. Also, with the in- 
crease of the spin of the black hole, it is seen that the 
disk-outflow coupled region shrinks considerably, attain- 
ing a steeper nature. Literally speaking, the spin of the 
black hole directly influences the nature and region of the 
outflow. The tendency of the outflow region to get con- 
tracted to the inner radius suggests that the outflow is 
more likely to exuberantly emanate from inner region of 
the disk for rapidly spinning black holes which have more 
gravitating power. Therefore, with the increase in spin, 
the disk-outflow region becomes more dense and more 
susceptible to eject the matter with a greater power. 

4. DYNAMICS AND NATURE OF THE FLOW 

As the system is gas pressure dominated and strongly 
advective, it is more receptive to strong outflows. Figures 
3 and 4 describe the variation of flow parameters as func- 
tions of radial coordinate r along the disk-outflow surface 
for various spin of the black hole. With the increase of 
the spin of the black hole, v r increases at the inner re- 
gion of the disk. It is seen that along the disk-outflow 
surface, v r becomes zero beyond a certain distance Rj S , 
which is also termed as zero v r surface illustrated in §3.1, 
which is explained it detail in the next section in order 
to understand the outflow power. The zero v r surface 
extends upto more inner region of the disk with the in- 
crease of spin of the black hole, indicating the fact that 
outflows and then jets originate from more inner region 
of the disk for rapidly spinning black holes. The increase 
of c s with spin (Fig. 4) indicates that the temperature 
of the disk-induced outflow is higher for rapidly spinning 
black holes. It is found from Fig. 4 that the maximum 
temperature of the accretion-induced outflow varies from 
~ 10 11 to 10 12 K corresponding to zero to maximal spin 
of the black hole. The truncation of the curves in the 
inner region symbolizes the inner boundary of the disk- 
outflow region as mentioned in §3.1. 
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Fig. 3. — Variation of radial velocity with radial coordi- 
nate. Solid, dashed, dot-dashed and dotted curves are for a = 
0,0.5,0.9,0.998 respectively. Other parameter 7 = 1.5. 
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Fig. 4. — Same as that of Fig. 3 but variation of sound speed. 
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5. ENERGETICS OF THE FLOW 

Accretion by a black hole is the primary source of en- 
ergy of the mass outflow from the inner disk region. The 
energetics of the accretion-induced outflow are mainly 
attributed to the mass outflow rate and the power ex- 
tracted by the outflow from the disk. The derivation of 
the mass outflow rate had been elaborated in Gf 0. We 
follow the same procedure here, however including the 
spin information of the black hole. The mass outflow 
rate is given by 



M 3 {r) 



Anrp(h surf )v z (h sur f )dr + . 



(19) 



where the constant Cj is determined by an appropriate 
boundary condition. From the adiabatic conditions, p 
can be written in terms of c s which is already deter- 
mined as a function of 7'. The corresponding proportion- 
ality constant is determined at a radius, outside which 
the contribution to the mass outflow rate is negligible 
(v z — 0). However, the total mass accretion rate M 
(which is sum of the inflow rate and the outflow rate) 
can be obtained by integrating the continuity equation 
along the radial and vertical directions. Hence, the con- 
stant is computed easily by supplying the values of v r 
and c s at that radius in the expression of M . As we have 
discarded any possible small amount of the outflow at 
that radius, the actual value of the constant, and hence 
the outflow power, could be slightly over estimated. 

Mj(r) in Eqn. (|f 9[) refers to the rate at which the 
vertical mass flux ejects from the disk-outflow surface. 
Figure 5 shows the variation of Mj profiles with the spin 
of the black hole. It is seen that with the increase of 
spin of the black hole, Mj increases. All the profiles 
have been shown considering a supermassive black hole 
of mass ~ I0 8 Mq with a mass accretion rate at infinity, 
M ~ 10~ 2 MEdd, where Mes,6, is the Eddington mass ac- 
cretion rate ~ 1.44 x 10 25 gm/s. Low mass accretion rate 
(sub Eddington accretion flow) is in conformity with our 
gas pressure dominated advective disk paradigm. It is 
found that Mj increases with the increase of mass of the 
black hole and mass accretion rate as well. The trunca- 
tion of the curves at an inner radius indicates the inner 
boundary of the disk-outflow coupled region, explained 
in detail in the next paragraph. 

In computing the power extracted by the outflow from 
the disk, we follow the same procedure as in G10. Thus 
the power of the outflow is given by 

p j( r ) = [ 4nr + — + <j>G)pv z dr, 

J IA 2 7- 1 P J J haurf 

(20) 

which is the total power removed from the disk by the 
outflow along the disk-outflow surface. In measuring the 
net power of any astrophysical jet, it appears that the 
computed power Pj will then be the initial power of the 
jet. Figure 6 depicts the variation of the power Pj with 
r for various spin parameters of the black hole. If we 
meticulously investigate the nature of the power profiles, 
we observe that when the radial distance is less than 
a certain value, Pj begins to decrease. The decreasing 
trend of Pj infers that the characteristic flow is bounded 
(integrand of Eqn. (|20p carrying the information of the 
vertical energy flux becomes negative). 




10 20 30 40 50 



Fig. 5. — Variation of mass outflow rate as a function of radial 
coordinate, when (a) a = 0, (b) a = 0.5, (c) a = 0.9, (d) a = 0.998. 
Other parameter 7 = 1.5. 
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Fig. 6. — Variation of outflow power as a function of radial 
coordinate, when (a) a = 0, (b) a = 0.5, (c) a = 0.9, (d) a = 0.998. 
Other parameter 7 = 1.5. 

This occurs owing to the fact that in the extreme in- 
ner region of disk, due to its strong gravitating power the 
starved black hole sucks all of the matter in its sphere 
of influence, even if there is any outflow emanating from 
the disk. We identify this inner transition radius as Rj t , 
beyond which no outflow occurs. We attribute Rjt as the 
inner boundary of the disk-outflow region. In all of the 
previous profiles, the said truncation of the curves is as- 
cribed to Rjt ■ We describe the dynamical variables in our 
study within the region between an outer boundary and 
Rjt , and within this prescribed region strong outflows are 
most plausible to originate, intrinsically coupled to the 
disk. The power profile shows that with the increase in a, 
Pj increases as well as Rjt gradually shifts to the vicinity 
of the black hole, indicating the fact that outflow region 
moves further inward. We show the variation of total Pj 
extracted from the disk-outflow region with the spin of 
the black hole in Fig. 7a. Considering a black hole of 
mass ~ I0 8 M Q , as seen in AGNs and quasars, accreting 
with M ~ l0~ 2 M E dd> it 

is seen that Pj ~ I0 41 erg js for 
a = 0. For a maximally spinning black hole (a = 0.998) 
with same parameters, the computed Pj ~ 10 43 erg/s. 
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Fig. 7. — Variations of (a) net outflow power in the units of 
10 43 erg s _1 , (b) inner transition radius of the outflow, (c) peak of 
disk-outflow surface, as functions of spin of the black hole. 

Thus, there is an increase of two orders of magnitude 
of Pj with the increase of spin of the black hole from 
to 0.998. In an earlier work, Donea & Biermann (1996) 
showed that the power extracted by the outflow /jet from 
the disk increases with a. However, they did not com- 
pute the power extracted from the disk explicitly. The 
numerical simulations by De Villiers et al. (2005), in 
a different accretion paradigm including magnetic field, 
also concluded that the jet efficiency is possible to in- 
crease with the increase of spin of the black hole from 
to 0.998. Figure 7b shows the variation of Rjt with a. In 
describing the disk-outflow surface in §3.1, we have ar- 
ticulated the impact of spin on the disk-outflow coupled 
region. The geometry of the surface for a particular a is 
identified with a parameter Rj S (see Fig. 2). We show 
exclusively the variation of Rj S with the spin of the black 
hole in Fig. 7c. which too reveals that the fast rotating 
black hole retracts the outflow region towards it; a sign 
of a pure relativistic gravitation. 

6. DISCUSSION 

Blandford-Znajek process (Blandford & Znajek 1977) 
is still one of the most promising mechanisms to drive 
powerful jets in AGNs and XRBs. Although the exact 
mechanism of formation of the jet in the vicinity of the 
black holes is still elusive, and whatever might be the 
reason for the origin of jet, the said work is significant 
mostly due to two underlying reasons: (1) extreme grav- 
ity is indispensable to effuse strong unbounded flows in 
the vertical direction from the inner region of the accre- 
tion disk, (2) the spin of the black holes, which is purely 
a relativistic effect, is essential to power strong outflows 
and jets. 

In the present study, we have neither laid importance 
to the nature of the outflow nor invoked any aspect for 
the origin of outflows or jets. Indeed, the distinctive or 
definitive understanding of the formation of strong out- 
flows or jets is till unknown. Notwithstanding, the most 
distinguishable and obvious picture in this case is that 
outflows and jets observed in AGNs and XRBs can only 
originate in an accretion powered system, where the ac- 
cretion of matter around relativistic gravitating objects 
like black holes acts as a source, and outflow and then jet 



takes the form of one of the possible sinks (the other sink 
is the central nucleus). The dynamics of the outflowing 
matter should then be intrinsically coupled to the accre- 
tion dynamics macroscopically through the fundamental 
laws of conservation (of matter, momentum and energy). 
The outflow is unbounded and the total energy just at 
the base of the outflow should be positive. The accretion 
flow should be bounded in the vicinity of the central ob- 
ject as the central potential is attractive. However, the 
strength of the unbounded flows in the form of jets is 
distinctly proportional to the attractiveness of the cen- 
tral gravitational potential field. This paradox is well 
manifested in the observed universe, as relativistic jets 
are more populated around extreme gravitating objects 
like black holes. Also noticeably, length scale of jets in- 
creases from microquasars to quasars, which are suppos- 
edly harboring stellar mass and supermassive black holes 
respectively. 

Thus in any theoretical modeling of the accretion and 
outflow, it can be presumably argued that the mathe- 
matical equations governing the dynamics of the inflow 
and outflow should inherently be correlated and evolved 
self-consistently without any ad hoc proposition, as ac- 
cretion and outflow should not be treated as dissimilar 
objects. Second, the relativistic gravitational effect of 
the black hole should be incorporated, as the nature of 
gravity is the cornerstone to both the accretion and the 
unbounded outflow. To capture this essential physics, in 
our present study to understand the connection between 
disk and outflow, we have incorporated the general rela- 
tivistic effect of the spinning black hole through a pseudo- 
Newtonian approach. Although pseudo-Newtonian for- 
malism is an approximate method to mimic the space- 
time geometry of the Kerr black hole, yet it captures the 
important salient features of the corresponding metric, 
and thus can be used to examine the nature of outflows 
from the inner region of the disk. 

In §2 we have described the general disk-outflow cou- 
pled hydrodynamic equations in the inviscid limit follow- 
ing GM09. The necessity to simplify our model equations 
for an inviscid flow is discussed in §2. Although GM09 
explored the 2.5-dimensional accretion-induced outflow 
for a fully viscous system, they used a self-similar ap- 
proach in order to solve the necessary partial coupled 
differential equations. In addition, they neglected the 
most indispensable effect of relativistic gravitation or, 
precisely, the effect of spin of the black hole. In the 
present paper, we have solved the disk-outflow model 
equations in a more general 2.5-dimensional paradigm, 
while trying to limit our assumptions to the least possi- 
ble extent theoretically. One of the important premises 
we have made is the relationship between v z and c s , 
which we have established empirically. We have not used 
the height integrated equations which are mostly valid in 
the circumstances where the dynamical fluid parameters 
are likely to be independent of z. Without presuming 
the fact that the outflow originates from the surface of 
the disk, as most of the authors do, we have logically 
constructed a disk-outflow surface with properly defined 
boundary conditions. One of the most important com- 
putations we have done is to evaluate the mass outflow 
rate and the power of the outflow extracted from the 
disk self-consistently in the inviscid limit, unlike the pre- 
vious works (e.g. Donea & Biermann 1996; Blandford 
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& Begelman 1999; Xie & Yuan 2008). We have found 
that the spin of the black hole plays a crucial role in de- 
termining the structure, dynamics and the energetics of 
the outflow coupled to the disk. With the increase of 
the spin, the outflow region extends further inward and 
then the disk-outflow region shrinks and compresses (see 
Fig. 2). As a result, the outflow and then any plausible 
jet is likely to eject out from an extreme inner region 
of the flow around a rapidly spinning black hole with a 
greater efficiency Note that the higher spin results in 
the system to get more compressed with a greater out- 
flow power. Therefore, the efficiency of outflow and jet 
is directly related to the disk scale-height and hence the 
disk-outflow surface. Previously, in a different context, 
Mckinney & Gammie (2004), while examining the elec- 
tromagnetic luminosity of a Kerr black hole, assumed the 
ratio of the disk height to the radius (h/r) constant, ir- 
respective of the black hole spin. However, as seen in the 
present work, the constant h/r for different spin of the 
black hole may not be an obvious choice. 

The power extracted by the outflow from the disk not 
only depends directly on the mass of the black hole and 
the initial mass accretion rate of the flow, but also on 
the spin of the black hole. With our model, keeping the 
black hole mass and the accretion rate the same, the 
power of the outflow increases with the spin of the black 
hole and the computed power differs in two orders of 
magnitude between non-rotating and maximally rotating 
black holes. We have restricted our study vertically up 
to the region where the inflow and the outflow are least 
coupled, i.e. the disk-outflow surface. Above this surface 
accretion ceased to exist and probably the outflow gets 
decoupled from the disk, accelerates and eventually forms 
relativistic jet. The modeling of the astrophysical jet is 
altogether a different issue and is beyond the scope of the 
present work. Nevertheless, it can be effectively argued 
that in modeling the dynamics of the jet, the computed 
outflow power may serve as an initial power fed to the 
jet. Thus, the dynamics and the energetics of the jet will 
eventually be related to the black hole spin. 

According to the unification scenario, it is possible 
to device a single astrophysical scenario, which can 
broadly explain the observed multitude of different types 
of AGNs (see e.g. Antonucci 1993; Urry & Padovani 
1995). Flat Spectrum Radio Quasars (FSRQs) and 
BL Lacs are probably the two most active types of 
AGNs which are collectively referred to blazar. It is ob- 
served that BL Lacs are relatively low luminous than FS- 
RQs. Although both of these galaxies show high energy 
emissions, their spectral properties are different (Bhat- 
tacharya et al. 2009) which indicate that they are differ- 
ent source classes. According to the unification scheme, 
for FSRQs the line-of-sight is almost parallel to the jet 
and, hence, strong relativistic Doppler beaming of the jet 
emission produces highly variable and continuum domi- 
nated emission. As one moves away from the jet axis, the 
central continuum emission falls and the nucleus looks 
like a FR-II galaxy. Similarly, BL Lacs are considered to 
be a sub class of FR-I galaxies whose line-of-sight is al- 
most parallel to the jet axis. The present work suggests 
that the total mechanical power of an outflow proportion- 
ately increases with the spin of the central supermassive 
black hole; higher the spin, stronger is the outflow. It is 
reasonable to consider that strong outflow can lead to a 



strong jet, and hence, one can expect to observe higher 
luminosity. Therefore, the work suggests that BL Lacs 
are slow rotators than FSRQs. 

For the theoretical formulation of the disk-outflow cou- 
pling even with approximations, one needs to be very 
thoughtful for the proper foundation of the model which 
essentially needs to solve hydrodynamic or magnetohy- 
drodynamic conservation equations in presence of strong 
gravity. The flow parameters vary in 2.5-dimension and 
are coupled to each other. Limited observational inputs 
put irremediable constraint on the boundary conditions 
as well as the fundamental scaling parameters, governing 
the coupled dynamics of the accretion and outflow. The 
inadequacy of an effective mathematical tool to handle 
partial coupled differential hydrodynamic equations for 
compressible flow motivates us to invoke approximations 
and assumptions. Despite of this fact, one needs to ex- 
plore the possibilities to examine the accretion-induced 
outflow. The question then arises: what are the valid as- 
sumptions and to which extent they can be considered? 
In the present work, we have addressed this question in 
the following way. 

(1) As the extreme gravity is the most important aspect 
to effuse jet from the disk, we have incorporated its effect 
through a pseudo-general-relativistic potential. 

(2) The disk and outflow should not be treated as dissim- 
ilar objects, and hence their correlated-dynamics should 
be essentially governed by the conservation laws. The 
energetics of the accretion-induced outflow would then 
be evaluated self-consistently as is shown in our work. 

(3) Any unbounded flow in the form of outflow is more 
plausible to emanate from a hot, puffed up region of 
the accretion flow (may be low/hard state of the black 
hole). Hence we have formulated our model in a 2.5- 
dimensional, strongly advective paradigm, surpassing the 
simplicity of height integration. 

(4) We have approximated our system to an inviscid 
limit, whose reason has been argued in §2. Nevertheless, 
in any future work, viscosity should be incorporated into 
the flow to make it more realistic. 

(5) We do not account for the mechanism of formation of 
strong jets, like due to magnetic field or strong radiation 
pressure, as the definitive mechanism of the formation 
of jet is still unknown. Indeed, it is beyond the scope 
to solve magnetohydrodynamic equations in the present 
scenario. However, assuming that the outflow resides, 
and is coupled to the disk, the governing conservation 
equations of matter, momentum and energy should be 
treated accordingly. 

(6) Power of the outflow and jet is expected to directly 
depend on the spin of the black hole. The spin, which is 
the signature of general relativity, can make an impact 
on the nature of the observed AGN classes, and thus the 
spin of the black has been incorporated in our study. As 
spin of the black hole has a direct impact on the flow pa- 
rameters, we obtain different outflow power for different 
spin. 

What can be the most defined parameter for the 
accretion powered system — mass accretion rate, mass 
of the central star/black hole or the spin of the central 
star/black hole? Our analysis has shown that the power 
extracted by the outflow from the disk proportionately 
increases with the spin of the black hole. It infers that 
if extreme gravity is essential to power the jet, then the 
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strength and the length-scale of the observed astrophys- 
ical jets directly depend on the spin of the black hole. 
If it is believed that the distant quasars harbor massive 
black holes of same mass scale and accrete matter with 
similar rate, perhaps, spin be the guiding parameter for 
the different observed AGN classes. We can end with 
the specific question: are the observed AGN classes 
astrophysical laboratories to measure the spin of the 
supcrmassive black holes? 

This work is supported by a project, Grant No. 
SR/S2HEP12/2007, funded by DST, India. The authors 
would like to thank the referee for making important 
comments which helped to prepare the final version of 
the paper. 
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APPENDIX 

Equation (|15[) consists of complicated terms, where 
A, B and C are given by following equations 
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